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Abstract 

We present a dynamical study of e^e^ and fi^n^ production in proton- 
nucleus and nucleus-nucleus collisions at CERN-SPS energies on the basis of the 
covariant transport approach HSD employing a momentum-dependent p-meson 
spectral function that includes the pion modifications in the nuclear medium 
as well as the polarization of the p-meson due to resonant p—N scattering. We 
find that the experimental data from the CERES and HELIOS-3 Collaborations 
can be described equally well as within the dropping p-mass scenario. Whereas 
corresponding dilepton q^^-spectra are found to be very similar, the inclusive 
dilepton yield in the invariant mass range 0.85 < M < 1.0 GeV should allow to 
disentangle the two scenarios experimentally. 
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1 Introduction 



The mechanism of chiral symmetry restoration at high baryon density is of fundamen- 
tal importance for the understanding of QCD [|T], 0, but a clear experimental evidence 
has not been achieved, so far. Enhanced antikaon yields at SIS energies - as compared 
to transport studies using bare kaon masses 0, ^ - point in this direction, but the 
present knowledge on the elementary production cross sections close to threshold - es- 
pecially for the AA^ production channel - does not yet allow for final conclusions. On 
the other hand, dileptons are particularly well suited for an investigation of the early 
phases of high-energy heavy-ion collisions because they can leave the reaction volume 
essentially undistorted by final-state interactions. Indeed, dileptons from heavy-ion 
collisions have been measured by the DLS collaboration at the BEVALAC [|, ^, |^ 
and by the CERES §, HELIOS [|, |ig, NA38 ^ and NA50 collaborations [g at 
SPS energies. 

The recent data on e+e^ or fi^fi^ spectra at SPS energies might provide valuable 
insight; Li et al. have proposed that the enhancement of e^e^ pairs in S + Au 
collisions as observed by the CERES collaboration ||^ stems from an enhanced p- 
meson production (via tt+tt" annihilation) due to a dropping p-mass in the medium. 
In fact, their analysis - which is based on an expanding fireball scenario - is supported 
by the microscopic transport calculations in ref. |T^. Meanwhile various authors 



have substantiated the observation of refs. ||T3|, |1J], i.e. that the spectral shape of 
the dilepton yield is incompatible with 'free' meson form factors [^, l7, 18, 19 . 
However, a more conventional approach including the change of the p-meson spectral 
function in the medium due to the coupling of p, vr, A and nucleon dynamics along the 
lines of refs. pO|, ^ |2^, ^ was found to be (roughly) compatible with the CERES data 
p3[] , as well, whereas the dimuon data of the HELIOS-3 collaboration ||^ could 



only be described satisfactorily when including 'dropping' meson masses |24, 25|. 
Meanwhile, our knowledge on the p spectral function has improved since - as first 
pointed out by Friman and Pirner - resonant p-A^ interactions in p-wave scattering 
significantly enhance the strength in the vector-isovector channel at low invariant 
mass. In fact, the CERES data for S + Au at 200 A GeV and Pb + Au at 160 A GeV, 
using an expanding fireball model, were found to be compatible within such a hadronic 
scenario P7[ . 

In the present work we report on novel developments in this direction by improving 
on both the reaction dynamics and the p spectral function {e.g. by constraints from 'jp 
data), which will be discussed in sect. 2. In sect. 3 the resulting dilepton spectra are 
then compared to the data from both CERES and HELIOS-3 collaborations simulta- 
neously; we also search for observables that might distinguish between the 'hadronic' 
and 'dropping mass' scenario, in particular dilepton g^-spectra and inclusive yields 
above the p mass. We finish with a summary and concluding remarks in sect. 4. 
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2 Heavy-Ion Reaction Dynamics and in-Medium 
Rho Propagator 

The dynamical description of proton-nucleus or nucleus-nucleus collisions is performed 
within the covariant transport approach HSD |^ which has been found to account 
for hadronic data from SIS to SPS energies reasonably well and has recently been used 
in the analysis of dilepton data W^, El, 12^. The dilepton production is calculated 



perturbatively by including the contributions from the Dalitz-decays rj 
uj — s> TT^l^l^, 7]' — s> ^l^l^ and the direct dilepton decays of the vector mesons p, 
and (f) where the p and the meson may as well be produced in vr+vr" and KK or 
Tip collisions, respectively. For a detailed description of all these processes we refer 
the reader to ref. |2^. We note that we discard baryon-baryon {BB), meson-baryon 
(mB) and meson-meson {mm) bremsstrahlung channels as well as the Dalitz decays 
of the baryon resonances in order to avoid double counting when employing the full 
p spectral function. 

The dilepton radiation from p mesons is calculated as 

dNi+i- 2M , , 

= — — im£'p(go, q; pb, T) (i) 

where Dp is the p-meson propagator in the hadronic medium, depending on the baryon 
density pB and temperature T as well as on energy go and 3-momentum g = |g| in 
the local rest frame of the baryon current ('comoving' frame). The invariant mass M 
is related to the p-meson 4-momentum in the nuclear medium as 

= ql - q' , (2) 

while 



1 
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is spin averaged over the longitudinal and transverse part of the p-propagator: 



ImD, = -(ImD^^ + 2ImDj) (3) 



T ni,T. rr.. ImS^-^(go, q; ps, T) 

ImiJ„' (qn,q:pB,T) = j-^ . 4 

The p-meson selfenergy S^'"^(go, <?; P_b, 7") is obtained by combining the effects of the 
different hadronic interactions as: 

= SpTTTT + ^pBB-^ + ^^ai + ^'^KK-t + ^p^Xi • (6) 

The explicit evaluation of the various selfenergy components is discussed in detail in 
ref. 12^. We here employ an improved version in the following respects: 



a) in addition to p-wave pN ^ B interactions {B=N, A, A^(1720) and A(1905)) 
we also account for s-wave excitations into iV(1520) and A(1700) resonances (as 
usual, the corresponding coupling constants are estimated from the pN partial 
decay width); 
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b) the coupling of the resonances to virtual photons is calculated within the im- 
proved vector dominance model of KroU, Lee and Zumino to avoid overes- 
timates of the B N'y branching ratios; 

c) the medium modifications of the two-pion selfenergy Sp^rTr are extended to arbi- 
trary 3- momentum within the recently developed model of Urban et al. . 

The combined p-meson selfenergy is then further constrained by experimental infor- 
mation on 7p and 7^4 absorption cross sections |^ which correspond to the limit 
0. The constraints by photo-production data have been stressed by several 



authors 



29, 19, 34 



35], ^] and thus qualify as stringent benchmarks for the p 
spectral function. We note that by resummation of hadronic interactions in the p- 
meson selfenergy Sp (eqs. ®,(@)) we incorporate dynamical effects which are not 
only linear in the baryon density pb-, but also of higher order. The resummation of 
interactions at fixed baryon density and temperature involves the assumption that all 
interaction terms add coherently - an assumption that might become questionable 



due to the short time scales of a few fm/c in heavy-ion collisions at SPS energies |l28 |. 
Our approach thus differs conceptually from the analysis performed in Refs. |36[| 



where only dilepton rates linear in have been considered and higher order coherent 
effects have been discarded. 

As an example - relevant for the heavy-ion reactions at SPS energies - we show 
in Fig. 1 the spin averaged —lmDp{q(),q;pB,T) as a function of the invariant mass 
M and the 3-momentum q for a temperature of 150 MeV at ps = 0,po,2po and 
3po, respectively. With increasing baryon density we find the p spectral function to 
increase substantially in width showing only minor structures at high density. Thus 
the lifetime of the p-meson in the nuclear medium becomes very short due to well 
established hadronic interactions. As a consequence, its average propagation in space 
is limited to distances less than 0.5 fm already at normal nuclear matter density (even 
for high momenta). 



The cross section for the pion annihilation channel tt 
in line with ref. 



TT 



e+e is taken 



37 as 



<J-n 



-(M) 



02 



[m 



bare\4 



) lmDp{qo,q;pB,T) 



(7) 



where k = (M^ — 4:m'lY^'^/2 is the pion 3-momentum in the center-of-mass frame and 
a is the fine structure constant. The pTrvr coupling constant and bare p mass are fixed 



to reproduce the pion electromagnetic form factor in free space 



In this way we ensure that the Fourier transformation of the isovector current- 
current correlation function in vacuum is consistent with the experimental information 
in terms of the vacuum p spectral function [35|. We note that constraints from chiral 



symmetry are not fully accounted for in our model for the spectral function since the 
ai(1260) meson is not treated as the chiral partner of the p meson |T9|. However, 
the dominant effects arise from p-wave irN- and resonant pA^ scattering, where chiral 
symmetry does not play a dominant role. 

Our dynamical calculations are carried out as follows: the time-evolution of the 
proton-nucleus or nucleus-nucleus collision is described within the covariant transport 
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approach HSD |2^, without dropping vector meson masses. Whenever a p-meson 
is produced in the course of the hadronic cascade (by baryon-baryon, meson-baryon 
or pion-pion colhsions), its 4- momentum in the local rest frame of the baryon current 
is recorded together with the information on the local baryon density, the local 'trans- 
verse' temperature and its production source. We note that the definition of a local 
temperature is model dependent; here we have used a logarithmic fit to the transverse 
Pt spectra of mesons at midrapidity. Without going into a detailed discussion of this 
issue we note that our results for dilepton spectra do not change within the numerical 
accuracy when using a constant temperature T = 150 MeV at SPS energies since 
ImDp depends rather weakly on T (at a given nucleon density). 



3 Comparison to HELIOS-3 and CERES Data 

We begin our discussion of the dilepton spectra at SPS energies with the p + W 
system at 200 GeV incident energy. Fig. 2 shows the spectral channel decomposi- 
tion as a function of the li^ ^~ invariant mass M in comparison to the data of the 
HELIOS-3 collaboration |^ including the experimental cuts in transverse mass and 
rapidity as well as the experimental mass resolution as in ref. [^. The spectrum 



is normalized to the charged particle multiplicity N^. in the pseudorapidity bin 3.7 

< 7] < 5.2. The upper part of Fig. 2 shows the dimuon spectra when employing the 
'free' p spectral function whereas the lower part displays our results when accounting 
for the full medium modifications of the p-meson propagator. As is apparent from 
Fig. 2, the spectrum for p -|- W can be fully accounted for by the electromagnetic 
decays of the r],rj' and vector mesons p^,uj and 0; contributions from meson-meson 
channels {tt^tt~ , KK,TTp) are of minor importance here. Furthermore, since the p- 
meson essentially hadronizes in the vacuum, there is no noticeable difference between 
the summed spectra in the upper and lower part of Fig. 2. 

The situation changes appreciably when turning to nucleus-nucleus collisions. In 
Fig. 3 we compare the results of our calculation for the differential dilepton spectra 
for S + W at 200 A- GeV with the experimental data ^ employing the same cuts 
and mass resolution as before. Again, the upper part shows our results for the 'free' 
p spectral function whereas the lower part represents our results including the full 
medium modifications of the p propagator. Even at forward rapidities as measured 
in central S -|- W (200 A-GeV) a sizable contribution for invariant masses 0.3 GeV 

< M < 0.7 GeV stems from the tt+tt^ ^ p — p+p^ channel. Also in the mass 
regime around 1 GeV there is a significant contribution from KK and vrp annihilation 
to dimuons which explains the enhanced 0/ (p + uj) ratio in S -|- W relative to p + W. 
As has been found earlier [^, ^ the experimental spectrum cannot be properly 
described in terms of a 'free' p spectral function: there is an excess for invariant masses 
M ^ 0.8 GeV and insufficient yield around M ^ 0.5 GeV. The description of the data 
is significantly improved when including the hadronic medium modifications of the p 
propagator (lower part of Fig. 3) - quite analogous to the 'dropping mass' scenario 



when employing a p-mass shift according to Hatsuda and Lee However, there 

persists a significant difference: whereas in the 'dropping mass' scenario (cf. Fig. 4 
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of ref. pj]) the dimuon yield is underestimated for M > 1.2 GeV the in- medium p 
spectral function seems to give sufficient pairs in this region due to its substantial 
broadening in dense matter (cf. Fig. 1). 

In Fig. 4 we also compare the results of our calculation for the differential dilepton 
spectra in Pb + Au at 160 A- GeV and 5 = 5 fm with preliminary experimental 
data |3^, The e^e'^ acceptance cuts in pseudorapidity (2.1 < rj < 2.65), transverse 
momentum {px > 0.175 GeV/c) as well as in the opening angle of the e~^e~ pair 
(0 > 35 mrad) are taken into account. Furthermore, the experimental mass resolution 
has been included in evaluating the final mass spectrum, which is normalized to the 
number of charged particles drich/dr] in the pseudorapidity bin 2.1 < i] < 3.1. Again, 
the upper part shows our results for the 'free' p spectral function whereas the lower 
part represents our results including the full medium modifications. For Pb + Au 
at 160 A- GeV (and semi central collisions) the dominant yield for invariant masses 
0.3 < M < 0.7 GeV stems from 7r"'"7r~ annihilation. In the mass regime around 
1 GeV there is again a large contribution from KK and vrp annihilation to dileptons. 
Within the present statistics, however, there is no unique signal for in-medium effects 
since the calculation with the 'free' spectral function (upper part) also describes the 
data except for one point at 0.7 GeV (cf. also ref. I]!^). The calculation with the 
full spectral function describes the data somewhat better to the same accuracy as 



within the 'dropping mass' scenario ||29|. We note that the comparison with the data 
in Fig. 4 has been performed for impact parameter 6 = 5 fm, for which our transport 
model gives a charged particle multiplicity dnch/drj ^ 260 in accordance with the 
experimental normalization. 

Since the present dilepton invariant mass spectra do not discriminate the different 
models we also explore if the momentum distributions of the lepton pairs give some 
further insight (as suggested by Friman and Pirner and in ref. [^). Our calculated 
transverse momentum spectra - integrated over rapidity and the invariant mass range 
0.3 < M < 0.7 GeV - are shown in Fig. 5 for Pb + Au at 160 A-GeV and 6 = 2 fm 
using gT^-bins of 50 MeV/c. The dashed line in Fig. 5 displays the result with the 
'free' spectral function while the solid line shows the result for the in-medium spectral 
function. Apart from an overall increase of the spectrum we do not observe significant 
changes of the slope in qt - as was the case for the 'dropping mass' scenario (cf. Fig. 13 
of ref. |2^). Thus, the qt spectra do not qualify in disentangling the different schemes 
either. 



Another signature studied in ref. |29] and also suggested in |36] is related to the 



total dilepton yield between the free u and peaks where the dominant contribution 
stems from p decays. In this regime the 'dropping mass' scenario leads to a reduction 
by about a factor 2-3 as compared to a 'free' p spectral function for Pb + Pb at 
160 A-GeV while a broadening of the p spectral function due to hadronic inter- 
actions shows an enhancement. This feature is demonstrated quantitatively in Fig. 6 
for Pb -|- Au at 6=5 fm for the different scenarios employing an experimental mass 
resolution of AM=10 MeV. We find a difference by about a factor of 5 between the 
'dropping mass' scenario (dashed line) and the hadronic spectral function approach 
(solid line) which might be accessible experimentally with the CERES detector at im- 
proved mass resolution. We note, furthermore, that in the spectral function approach 
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the relative enhancement is most pronounced at M ~ 0.3 GeV while in the dropping 
mass scheme we find a maximum enhancement at M ~ 0.5 GeV for this system. 



4 Summary 



On the basis of the covariant transport approach HSD ^] we have studied dilep- 
ton production in proton and heavy-ion induced reactions at SPS energies. Vari- 
ous contributions are accounted for: the Dalitz-decays of rj.oj.rj' mesons, the direct 
dilepton decays of the vector mesons p, uj and (f) as well as KK and vrp channels. 
It is found that for p + W collisions at 200 GeV the mesonic decays almost com- 
pletely determine the dilepton yield, whereas in S + W and Pb + Au reactions the 
7r"'"7r~ — >■ p — >■ KK — > — >• /+/~ annihilation channels and vrp collisions con- 

tribute substantially. The experimental data taken by the HELIOS-3 and CERES 
collaborations P, |4y] are in general underestimated by the calculations for invari- 
ant masses 0.35 GeV < M < 0.65 GeV when using a 'free' p spectral function (in 
accordance with refs. 



13, 15, E 



It turns out that both the 'dropping mass' scenario - using either Brown/Rho 



scaling |T^, |25| or dropping vector meson masses according to QCD sum rule studies 
[ll4| , |4[ ^ - as well as the hadronic spectral function approach lead to dilepton 
spectra that are in good agreement with the experimental data for all systems at SPS 
energies (cf. also ref. [^). In order to discriminate between the different models we 
have investigated dilepton transverse momentum distributions for invariant masses 
0.3 < M < 0.7 GeV in Pb -|- Au collisions; unfortunately, only little differences have 
been found as compared to the dropping mass scenario . However, the in-medium 
broadening of the p spectral function results in a much higher dilepton rate than 
the dropping mass scenario above the tu-mass. Therefore one might experimentally 
distinguish the two different scenarios by the total dilepton yield in the invariant 
mass range 0.8 < M < 1 GeV in Pb -|- Au collisions (cf. Fig. 6), provided that the 
experimental mass resolution is on the order of about 10 MeV. 
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Figure Captions 



Fig. 1: The (negative) imaginary part of the p propagator averaged over the longitudi- 
nal and transverse components as a function of the invariant mass M and the 
momentum q for baryon densities of 0, 1 po, 2 po, and 3 po and temperature 
T = 150 MeV. Note the different absolute scales in the individual figures. 

Fig. 2: The calculated dimuon spectra (full solid line) for p + W at 200 A-GeV in 
comparison with the data from ref. [^. The thin lines indicate the individual 
contributions from the different production channels including the HELIOS- 
3 acceptance and mass resolution; i.e. starting from low M: rj —>■ 7e"*"e~ 
(dashed line), u — > 7r°e+e~ (dot-dashed line), r]' •ye'^e' (long dashed line); 
for M ^ 0.8 GeV: u e~^e~ (dot-dot-dashed line), p° e+e^ (dashed line), 
7r+7r~ — > p — > e+e^ (dot-dashed line); for M 1 GeV: e+e~ (dot-dashed 
line), vrp e~^e~ (dot-long dashed line), KK e"^e~ (dashed line). The upper 
part shows the result for a 'free' p spectral function whereas the lower part is 
obtained for the full p spectral function. 

Fig. 3: Dimuon invariant mass spectra for central collisions of S + W at 200 A-GeV 
(full solid lines) in comparison to the data of the HELIOS-3 collaboration 0. 
The upper part shows the results of a calculation with a 'free' p spectral function 
whereas the lower part includes the full p spectral function as described in the 
text. The assignment of the individual contributions is the same as in Fig. 2. 

Fig. 4: Dielectron invariant mass spectra for semi central collisions of Pb + Au at 
160 A-GeV (full solid lines) in comparison to the preliminary data of the CERES 
collaboration |^. The upper part shows the results of a calculation with a 'free' 
p spectral function whereas the lower part includes the full p spectral function 
as described in the text. The assignment of the individual contributions is the 
same as in Fig. 2. 

Fig. 5: The dielectron transverse momentum distribution for the invariant mass range 
0.3 < M < 0.7 GeV for Pb + Au at 160 A-GeV at 6 = 5 fm for the 'free' p 
spectral function (dotted line) and with our in-medium spectral function (solid 
line) . 

Fig. 6: The dielectron invariant mass spectra for Pb + Au at 160 A-GeV at 6 = 5 fm 
for the 'free' p spectral function (dotted line) and with our in-medium spec- 
tral function (solid line). The dashed line corresponds to the 'dropping mass' 
scenario investigated in ref. PP|. 
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Fig. 1 
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Fig. 5 
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